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ABSTRACT 



We have observed the massive star forming region associated with the early 
B protostar G192.16-3.84 in NH3(1,1), 22.2 GHz H2O masers, 1.3 cm continuum 
emission, and at 850 /xm. The dense gas associated with G192.16 is clumpy, 
^ I optically thin, and has a mass of 0.9 M©. The NH3 core is gravitationally 

00 I unstable which may signal that the outflow phase of this system is coming 

to an end. Water masers trace an ionized jet 0.8" (1600 AU at a distance of 
O ; 2kpc) north of G192.16. Masers are also located within 500 AU of G192.16, 

^ I their velocity distribution is consistent with but does not strongly support the 

(--| ! interpretation that the maser emission arises in a 1000 AU rotating disk centered 

0-(\ on G192.16. Roughly 30" south of G192.16 (0.3 pc) is a compact, optically thick 

O I ('?" = 1-2) NH3 core (called G192S3) with an estimated mass of 2.6 Mq. Based 

'^ I on the presence of 850 /^m and 1.2 mm continuum emission, G192S3 probably 

^ ' harbors a very young, low-mass protostar or proto-cluster. The dense gas in the 

G192 S3 core is likely to be gravitationally bound and may represent the next 
rS ' site of star formation in this region. 
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1. Introduction 



The wealth of observational evidence suggests that outflows from low-mass protostars 
are almost certainly powered by disk accretion where magnetic stresses control the balance 
between inflow and outflow and determine the outflow structure far from the protostar. 
Intermediate mass protostars, up to early B stars, can probably also be described by this 
general picture however the details are not as clear. The presence of accretion disks and 
outflows traced to within 50-100 AU of early B protostars suggests that they are also likely 
formed via accretion (e.g.. Shepherd 2003, Chini et al. 2004). 

The structure and kinematics of circumstellar material around early B protostars 
directly affects the formation and evolution of the central object (s). Circumstellar disks 
provide a reservoir of material that can transport angular momentum away from the central 
object in two possible ways: 1) disks create outflows which transport angular momentum 
(e.g., Shu et al. 2000; Ouyed & Pudritz 1999); and 2) if M^isk/M^ > 0.3 the disk may be 
locally unstable and provide an additional means of transporting angular momentum away 
from the protostar (e.g., Laughlin & Bodenheimer 1994; Yorke, Bodenheimer, & Laughlin 
1995; Shepherd, Claussen, & Kurtz 2001, hereafter SCKOl). 

Massive star formation tends to occur in compact clusters and it can be difficult to 
separate individual outflows from each other (e.g., Beuther, Schilke, & Stanke 2003; Beuther 
& Schilke 2004). Even when a single, early B protostar dominates the flow dynamics in a 
cluster, companions or close members of a cluster can interact with an accretion disk causing 
it to precess or be truncated (e.g., Fendt & Zinnecker 1998 and references therein; Shepherd 
et al. 2000). The strong radiation field from massive protostars and cluster members can 
also significantly affect the cluster and disk evolution (e.g., Konigl 1999; Yorke & Sonnhalter 
2002). 

Understanding the interactions between cluster members, dense gas, and the radiation 
field is important to be able to unravel the processes that affect the formation and evolution 
of massive star forming regions. One early B protostar that may shed some light on this 
issue is G192. 16-3.84 (hereafter G192.16, also known as IRAS 05553+1631). G192.16 is the 
most massive object in a protostellar cluster at a distance of 2 kpc. An ultracompact (UC) 
Hii region with luminosity ~ 3 x 10^ L© implies the presence of an early B star with a mass 
of 8 to 10 M© (Hughes & MacLeod 1993; SCKOl). The detection of NH3(1,1) (Molinari 
et al. 1996), strong 22 GHz H2O maser emission (Codella, Felli, & Natale 1996), and 1 
& 3 mm continuum emission (Shepherd et al. 1998; Beuther et al. 2002) also suggests 
that this is a site of recent star formation. G192.16 is at the center of a 100 Mq molecular 
outflow with mass outflow rate Mfiow ~ 10""^ Mq yr^-*^ and dynamical age ~ 2 x 10^ years 
(Snell, Dickman, & Huang 1990; Shepherd et al. 1998, Devine et al. 1999). Recent 2.1 /im 
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observations detect a very red object near the position of the radio-detected UCHii region. 
If the 2.1 /im source is a direct detection of G192.16, then the protostellar photosphere is 
extincted by Ay = 29 — 33 magnitudes (Indebetouw et ah 2003). 

Based on 22.2 GHz H2O maser and centimeter continuum observations, Shepherd & 
Kurtz (1999) proposed that the G192.16 massive young stellar object (YSO) is surrounded 
by a 1000 AU rotating, molecular torus. The size of this molecular torus or possible disk is 
smaller than the handful of structures detected around other early B protostars. For example, 
IRAS 20126+4104 appears to have a 1700 AU flattened, rotating torus detected in CH3CN 
(Cesaroni et al. 1999) and a similar sized structure is detected in NH3 emission around 
AFGL 5142 (Zhang et al. 2002). It is unclear whether the circumstellar material detected in 
these dense core tracers are flattened, rotating cores, tori, or large disks. Significantly larger 
circumstellar structures are also observed: a protostar in M17 appears to have a torus or 
fiared disk of diameter 20,000 AU as mapped with near-infrared imaging (Chini et al. 2004); 
and fiattened, rotating molecular cores are inferred from high density tracers such as H^'^CN 
(NGC 7538S; Sandell, Wright, & Forster 2003) and NH3 (IRAS 20126+4104; Zhang, Hunter 

6 Sridharan 1998). 

Approximately 1" north of G192.16, H2O masers are also detected along a strip of 
ionized gas with a negative spectral index implying the presence of a mixed synchrotron and 
thermal emission jet (Shepherd & Kurtz 1999). The driving source for this collimated jet 
has not been detected. 

Continuum observations at 7 mm by SCKOl with a resolution of ~ 40 mas did 
not detect the 1000 AU torus predicted by Shepherd & Kurtz (1999). Instead, compact 
warm dust emission was detected within 50 mas of the massive YSO. A model fit to the 

7 mm continuum emission suggested an 8 M© central protostar (G192S1) with a possible 
protostellar companion (G192S2), a 130 AU diameter disk with a mass of 3-20 M©, and an 
outfiow with ~ 50 degree opening angle. The non-detection of emission from the the 1000 AU 
torus by SCKOl could be the result of resolving out the torus with the 7mm observations 
or the dust temperature could be too low to emit significantly in 7 mm continuum. Recent 
observations at 450 and 850 /im (at resolution 8" and 14", respectively) show that there is 
extended warm dust emission centered on G 192. 16 (Williams, Fuller, & Sridharan 2004). 
The far-infrared emission is centrally peaked and corresponds to a gas+dust mass of roughly 
26 to 64M0 (Williams et al. assumed a kinematic distance of 2.5 kpc; at a distance of 2kpc 
the mass would be 16 to 4OM0). 

In an effort to understand the cluster and circumstellar environment of G192.16, we 
present observations of NH3(1,1) tracing dense gas and H2O masers near the massive 
protostar, along with archive data of 850 /im continuum emission. 
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2. Observations 

2.1. Very Large Array 

Observations of NH3(1,1) line emission at 23.694 GHz were made with the National 
Radio Astronomy Observatory's"^ Very Large Array (VLA) in "D" configuration. Baselines 
between 35 m and 1.03 km could detect a largest angular emission scale, Olas ~ 60". The 
estimated uncertainty of the flux calibration is ~ 3%. Calibration and imaging was performed 
using the AIPS++ data reduction package. The continuum emission was subtracted in the uv 
plane. The line-only emission was imaged using Briggs robust uv weighting and deconvolved 
with CLEAN. Continuum-only emission was recovered and imaged separately using natural 
uv weighting and CLEAN deconvolution. 

Observations of H2O masers at 22.235 GHz were made with the VLA in "B" conflgura- 
tion. Calibration and imaging was performed using the AIPS data reduction package. The 
data were imaged using natural uv weighting and CLEAN deconvolution. A summary of the 
VLA observations is given in Table 1. 

Table 1: VLA Observation Summary 





NH3(1,1) 


1.3 cm continuum' 


H2O maser 


Observation Date 


2003 April 25 




2001 June 17 


Configuration 


D array 




B array 


Rest Frequency (GHz) 


23.6945 




22.2351 


Beam Size (arcsec) 


2.93 X 2.55 


2.92 X 2.55 


0.68 X 0.27 


Beam P. A. (degrees) 


44.8 


44.8 


-61.65 


Resolution ( km s^-*^) 


0.31 




0.33 


Total Bandwidth ( km s^^) 


39.4 




84.3 


Gain Calibrator 


J0530+135 




J0530+135 


Flux Calibrator 


J0137+331 (3C48) 




J0137+331 (3C48) 


Number of Channels 


128 


1 


256 


Peak Flux ( mJy beam^^) 


25.3 


2.8 


7,460. 


Map RMS ( mJy beam"^) 


2.2 


0.39 


28 



T Continuum emission subtracted from NH3 line in the uv plane. Observational parameters are the 
same as for NH3. 



''The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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2.2. Very Long Baseline Array 

Observations of Ii20 masers were made at tliree epoclis witli a 2-4 montli time 
separation using the Very Long Baseline Array (VLBA). Observations were not phase- 
referenced, thus absolute position information was lost. Calibration and imaging was 
performed using the AIPS data reduction package. The data were imaged using natural 
uv weighting and the image was deconvolved with CLEAN. Masers were detected in each 
observation and were quite variable with flux densities ranging from 15 mJy beam~^ to 
230 Jy beam^^. A summary of the VLBA observations is given in Table 2. 

In all three VLBA epochs, masers were present at two locations in the image. The epochs 
were aligned assuming the brightest maser spot remained constant in position. Based on 
the scatter distribution of maser features, the relative positions between epochs are probably 
accurate to within 5 mas. Because absolute position information was not available, proper 
motions of individual masers could not be determined. 

An estimate of the absolute position of the VLBA masers was obtained by comparing 
the brightest VLBA H2O masers observed on 2001 May 01 with the brightest VLA masers 
observed on 2001 June 17 at similar velocities. Although the VLA and VLBA observations 
were not obtained simultaneously, the positions of three masers with relative positions 
(270 mas, -495 mas), (290 mas, 175 mas), and (Omas, Omas) remained roughly constant. 
The three spots from each observation were aligned to within 5 mas. Thus, the relative error 
between the VLA and all VLBA images is 7 mas (adding positional errors in quadrature) 
while the absolute error of the VLA image is ~ 0.1". 



Table 2: VLBA Observation 


Summary 








Epoch 1 


Epoch 2 


Epoch 3 


Observation Date 


2001 May 04 


2001 July 02 


2001 November 17 


Rest Frequency (GHz) 


22.2351 


22.2351 


22.2351 


Beam Size (mas) 


0.86 X 0.36 


1.10 X 0.33 


0.96 X 0.36 


Beam P. A. (degrees) 


-5.30 


-15.35 


-13.09 


Resolution ( km s""*^) 


0.21 


0.21 


0.21 


Total Bandwidth ( km s^^) 


53.9 


53.9 


53.9 


Bandpass Calibrator 


J0555+3948 


J0555+3948 


J0555+3948 


Number of Channels 


256 


256 


256 


Peak Flux ( Jy beam~^) 


8.7 


0.5 


228.7 


Map RMS ( Jy beam^^) 


0.007 


0.002 


0.002 
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2.3. James Clerk Mcixwell Telescope 

Submillimeter continuum observations were obtained from the Canadian Astronomy Data 
Centre Archive of SCUBA (Submilhmeter Common User Bolometer Array, Holland et al. 1999). 
Observations from Williams et al. (2004) on the night of 2000 March 18, were combined with earlier 
observations from the night of 1997 October 18 to increase the sensitivity in the final image. Three 
850 //m jiggle maps were reduced using the SURF (SCUBA User Reduction Facility) reduction 
package (Jenness & Lightfoot 1998). Each map was flat-fielded and corrected for chop throw (120"), 
extinction, and sky noise. The extinction was determined by using the r225 opacities determined at 
225 GHz with the tipper located at the Caltech Submillimeter Observatory and scaling this value 
to the 850 /im optical depth, rgso, using the relationship derived by Archibald et al. (2002). The 
final 850 ^um image was made by adding together the individual images with pixels sampled at the 
Nyquist limit (7"). The final image has a resolution of 15". The image was calibrated using Uranus 
jiggle maps observed within a few nights of G192.16. The total estimated flux uncertainty is ~30% 
based on comparisons of the peak flux in each separately calibrated image. 



3. Results 

Figures 1 &: 2 present an overview of the interferometric observations of NH3(1,1), 1.3cm 
continuum emission, and H2O masers observed with the VLA and VLBA. G192.16, detected in 
1.3cm continuum emission at 05*^58^13.53**, 16°3l'58.3" (J2000), is located within a region of 
clumpy NH3 emission. H2O masers are only detected within 1" of the massive protostar. A 
compact NH3 core, G192 S3^ is discovered ~ 30" south of G192.16 at 05'^58™ 13.67", 16°31'31.8" 
(J2000). No centimeter continuum emission or H2O masers are detected near G192S3. An image 
at 850 /im shows a strong central peak at the position of G192.16 along with an extension toward 
G192S3 (Fig. 3). 



3.1. Dense molecular gas and embedded protostars 

NH3 emission near G192.16 is clumpy and covers a region of roughly 20" x 10" (~ 0.2 x 0.1 pc 
at a distance of 2kpc). Figure 1 shows that the G192.16 protostar is located near the center 
of the clumps. The emission is detected in each channel at a level of only 2 — 3cj (Fig. 4). 
Spectra averaged over G192.16 and G192S3 (Fig. 5) show that the G192.16 NH3 emission peaks 
at vlsr = 5.7 km s^"*^ while G192 S3 is shifted by 1 km s^^ to 6.7 km s^^. The spectra show both 
the main inversion line and the first two hyperfine components of the NH3(1,1) lines, although the 



^SCKOl designated G192 SI as the B2 protostar and G192 S2 as the proposed companion located ^ 80 AU 
north of G192S1. Here, we use the term G192.16 to refer to the S1+S2 system. 
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satellite lines toward G192.16 are detected at a level of only Sa. Gaussian fits to the main lines 
result in a peak of 4 mJy beam"^ and FWHM of 1.52 km s^^ for G192.16 and 13 niJy beam^^ and 
FWHM of 0.77 km s~^ for G192 S3. Figure 6 shows a first moment map of the intensity-weighted 
velocity field of the NH3 emission, [/ I{v)v dv] / [J I{v)dv]. There is no obvious velocity gradient in 
either region although the narrow line width could easily prevent detection of a velocity gradient 
inG192S3. 

The relative intensities of the main and satellite components of the (1,1) transition can be 
used to estimate the optical depth of the main component, r, from the equation (Ho & Townes 
1983): 

ATX{J,K,m) _ 1 - e-^(-^'^'"^) 

ATX{J,K,s) ~ 1 - e-'^^iJ^K'"^) ^^' 

where T^iJ, K, m) is the measured antenna temperature of the (J,K) transition for the main line, s 
refers to the satellite line, and a is the satellite to main intensity (0.28 & 0.22 for the (1,1) transition). 
The ratio of intensity of the (1,1) and (2,2) transitions is sensitive to level populations and thus acts 
as a temperature probe for cold molecular gas. Unfortunately, we were not able to observe the (2,2) 
transition. Temperatures of cores in Orion range from 15 to 40 K with 30 K being typical of warm 
cores in the Orion NH3 filaments (Wiseman & Ho 1998). Thus, we assume a temperature of 30 K. 
Following the procedures and assumptions described in Ho & Townes (1983) and Wiseman & Ho 
(1998) and assuming an [NH3]/[H2] abundance ratio of 3 x 10~® (Harju, Walmsley, & Wouterloot 
1993), we estimate the NH3 column density, H2 number density, and total gas mass of the NH3 
cores (Table 3). The last column in Table 3 gives the range of values found in NH3 cores in Orion. 
Within the large uncertainties of our assumptions, the clump properties in the G192.16 star forming 
region are similar to those found in Orion. 

Table 3: Properties of the ammonia cores 









Comparison with 




G192.16 


G192S3 


NH3 Cores in Orion''" 


r (NH3 main component) 


0.06 


1.2 


0.1-2.1 


NH3 Column Density (cm^^) 


8.1 X 10^3 


8.2 X 10^4 


1.9 X 10^2 _ g_3 X ioi4 


H2 Number Density (cm^^) 


8.3 X 10^ 


1.6 X 10^ 


5.0 X 10^ - 6.3 X 10^ 


Total Gas Mass (M©) 


0.9 


2.6 


0.1-1.9 



t Wiseman k Ho (1998) 



Ionized gas is detected in free-free continuum emission at 1.3 cm at the position of 
G192.16. The source is unresolved with a flux density of 2.8 mJy (7(t). The number of Lyman 
continuum photons required to produce the central peak is log A'^i^ = 44.87 photons s^^, which 
corresponds to a single zero-age-main-sequence B2 star with Li^oi ~ 2.8 x 10^ Lq (Thompson 
1984). This is consistent with the luminosity derived by Shepherd et al. (1998) and Shepherd 
& Kurtz (1999) to within the uncertainties. A spectral energy distribution (SED) of G192.16 



(Fig. 7) shows our data in relation to previous measurements. The fit to the free-free 
emission at centimeter wavelengths has a slope of 0.3 ± 0.07 (i.e., Sy oc z/*^'^), consistent with 
a moderately optically thick Hii region. No free-free emission was detected toward G192S3. 



3.2. Water maser emission 

VLA observations of G 192. 16 show a field rich in H2O maser emission. Red-shifted 
masers located 0.1" to 0.2" (2000 - 4000 AU) E-NE of G192.16 have velocities between 6 
and 9 km s~^ relative to vlsr- Red and blue-shifted masers are also located along a stripe 
to the north of G192.16 where Shepherd & Kurtz (1999) discovered an ionized jet of mixed 
synchrotron and thermal emission. The masers associated with this jet tend to have red- 
shifted velocities to the NE and blue-shifted velocities to the SW along the jet axis. Figure 
1 shows the locations and velocities of the H2O maser emission detected in our observations 
as well as those detected previously. 

VLBA H2O maser images were aligned with the VLA maser images as described in 
Section 2.2. The maser emission is variable in position, intensity, and velocity although 
masers exist in every epoch. VLBA H2O masers are detected near the end-points of the 
synchrotron-thermal jet in each of the three epochs (Fig. 2). Masers associated with the jet 
tend to be red-shifted to the NE and blue-shifted to the SW and have similar velocities to 
those observed with the VLA. During the first epoch of VLBA observations, one position 
with red and blue-shifted maser emission was detected 100 mas (200 AU) NE of G192.16. 

Although every VLA observation (this work and previous observations) detected H2O 
masers near the massive protostar, two of three VLBA epochs did not detect emission in this 
region. H2O masers are notoriously variable and those in G192.16 are no exception. Thus, 
it could be that the region near G192.16 did not have masers during our last two VLBA 
observations. It seems unlikely that the VLBA would resolve out the emission given that 
the ~ 1 mas resolution of the VLBA observations corresponds to a size-scale of about 2 AU 
at a distance of 2kpc which is larger than typical H2O maser features found toward other 
massive star formation regions. However, we cannot rule out this possibility. 



3.3. Imaging at 850 /im 

The 850 /im image (Fig. 3) has a strong peak (S/N ~ 96) centered on the G192.16 
UGH II region (shown as a cross in Fig. 3) along with an extension to the south where 
G192S3 is detected in NH3 emission (shown as a plus symbol in Fig. 3). This extension is 



-9- 



visible in each of the three individual jiggle maps and is thus likely to be real. A similar 
extension of the emission from G192.16 to G192S3 can be seen in the 1.3 mm continuum 
image of Beuther et al. (2002) and in the 21.3 /im image taken by the Midcourse Space 
Experiment (MSX, figure not shown) although both data show this feature only at the 
2 — 3(7 level. The 10% contour level in Fig. 3 extends nearly an arcminute south of G192.16 
suggesting the presence of molecular gas up to 30" beyond G192S3. 

The mass of G192 S3 may be determined from the flux in a small aperture (15") centered 
on the NH3 peak position; however, this flux will be contaminated by the extended emission 
associated with the massive molecular core surrounding G192.16. To minimize this affect, we 
constructed a model image of the G192.16 core from an azimuthally averaged radial intensity 
profile. In creating this average intensity profile a 60° wide sector centered on G192 S3 was 
omitted to ensure that emission from G192 S3 was not included in the model. The resulting 
model was subtracted from the 850 /xm image to produce the image shown in the right panel 
of Fig. 3. The flux density of G192 S3 in a 15" aperture measured in this model-subtracted 
image is 130 ± 45 mJy beam"^. Assuming: (1) a distance of 2 kpc; (2) a 1:100 dust to gas 
mass ratio; (3) a total gas opacity at 850 ;um of 0.034 cm^g~^ appropriate for coagulated 
dust grains with thin ice mantles (Table 1, column 5 of Ossenkopf & Henning 1994); and (4) 
an average dust temperature of 30K; then the total gas+dust mass of the southern core is 
1.3M0. This mass is a factor of 2 lower than the mass determined from NH3 observations. 
The Ossenkopf & Henning opacities provide a good flt to the observed SEDs toward star- 
forming cores (e.g., Evans et al. 2001; Shirley, Evans, & Rawlings 2002; Young et al. 2003) 
but are probably uncertain to within a factor of a few. Also, the systematic uncertainty in 
aperture photometry performed on the subtracted model is probably ~ 50% . Given these 
uncertainties, the mass estimates based on NH3 and 850 /xm imaging agree remarkably well 
and suggests that they are sampling similar material. 

The peak flux density of the G192.16 core is 2.10±0.63 Jy beam~^ and the integrated flux 
density is 7.4 ±2.3 Jy in a 120" aperture (uncertainties include contributions from systematic 
flux calibration). The derived flux density agrees to within the errors with Williams et al. 
(2004). The Williams et al. flux densities at 450 and 850 /im and our 850 /xm point are 
shown in the SED in Fig 7. Using the same assumptions given above, the total mass is 
75M0. This is signiflcantly larger than the mass estimate based on NH3 emission (O.9M0) 
and is consistent with the interpretation that the NH3 emission results from only a small 
fraction of molecular material that remains in dense, compact structures. In contrast, the 
850 ;um mass estimate of 75M0 agrees to within the errors with the II7M0 derived from 
Ci8o(J=l-0) (Shepherd & Kurtz 1999). 
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4. Discussion 

4.1. Water maser activity 

It is well known that H2O masers are commonly observed in both low- and high-mass 
star forming regions and the presence of H2O masers has been linked to the evolutionary 
status of star formation in a region (e.g., Tofani et al. 1995, Kurtz et al. 2000, Codella et 
al. 2004). Most often, masers are associated with jets and outflows (e.g., Claussen et al. 
1998; Torelles et al. 2001, 2003; Moscadelli, Cesaroni, & Rioja 2000; Furuya et al. 2003, to 
name a few) although there is some evidence that H2O masers may trace circumstellar disks 
around massive protostars (e.g., Cesaroni 1990; Torrelles et al. 1997; Molinari et al. 2000; 
Seth et al. 2002; Goddi et al. 2004). 

Water masers north of G192.16 appear to trace a collimated jet although the actual 
driving source of the jet has not been detected. The velocity distribution of masers along 
the jet axis (e.g., mostly red-shifted to the NE and mostly blue-shifted to the SW) suggests 
that the driving object lies near the center of the masers. 

Shepherd & Kurtz (1999) proposed that the H2O masers near G192.16 arose in a 
1000 AU flattened, rotating disk surrounding the protostar. The velocities were consistent 
with Keplerian rotation and with the rotation axis of a C^^O core. Their proposed 1000 AU 
disk is shown as a grey torus in Figs. 1 & 2. Our new observations of H2O masers with the 
VLA and VLBA are consistent with this interpretation however they do not shed a signiflcant 
amount of light on the issue given that the VLBA observations only detected a single maser 
spot with both red and blue-shifted emission peaks within the proposed disk and the VLA 
observations only detected a hand-full of red-shifted masers E-NE of 0192.16. Thus, we can 
neither conflrm nor deny the proposal that a 1000 AU disk surrounds the 0192.16 protostar. 



4.2. Molecular core stability 

Two ammonia cores with very different characteristics are detected toward the 0192.16 
star forming region. The core surrounding the 0192.16 massive protostar is clumpy, roughly 
O.lSpc in diameter (assuming a distance of 2kpc), and has a mass of about O.9M0. The 
NH3 line is optically thin with a full width at half maximum (FWHM) of 1.5 km s~^ and 
the velocity fleld is somewhat chaotic. The mass estimate based on 850 /xm emission from 
warm dust is 75M0 and compares reasonably well with estimates based on O^^O observations 
(II7M0). The large discrepancy between the mass derived from NH3 and masses derived 
from 850 yum and C^^O suggests that more than 99% of the mass of the cloud resides in 
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material with a density less than ~ 10^ cm 

In contrast, G192S3 is compact (diameter = 0.05 pc), optically thick, and has a mass 
of 2.6 Mq. The line FWHM is roughly half that of the G192.16 core (0.77 km s'^); no 
velocity gradient is seen across the core although a gradient could easily be missed with the 
0.31 km s~^ resolution of our observations. The similar mass estimates from NH3 (2.6M0) 
and 850 /im emission (l.SM©) suggests that both methods may probe the same region. 

The gravitational binding energy of a cloud core is given by GM^i^^^/ciRdoud where 
Mcioud is the mass of the cloud, Rdoud is the radius, and ci is a constant that depends on 
the mass distribution (ci = 1 for p ex r^^). This can be compared to the cloud turbulent 
energy to determine if the cloud is likely to collapse or disperse. Following Myers (1983), 
the one-dimensional turbulent velocity dispersion a is related to the line width Af by 

8/n(2) m ^ ' 

where m is the NH3 molecule mass {17 aranxiriH where rriH is the mass of the hydrogen 
atom) and we take T = 30K. The cloud turbulent energy is then ^Mdoud o"^- A comparison 
of the gravitational binding energy and the turbulent energy of the G192.16 and G192 S3 
cores is given in Table 4. The ratio of Eturb/Egrav is 11.7 for G192.16 and 0.3 for G192 S3 
suggesting that G 192 S3 is gravitationally bound and will eventually collapse to form a 
protostar while the NH3 gas surrounding G192.16 is gravitationally unstable and will likely 
disperse. Although there are numerous uncertainties associated with these estimates (see 
previous section), the values are consistent with the scenario that the G192.16 protostellar 
outflow and ionizing wind may have significantly disrupted the natal molecular cloud and 
are in the process of dispersing the remaining dense molecular gas. G192S3, 30" (0.3 pc) 
south of G192.16, appears to be far enough away from the disrupting influence of the massive 
protostar to allow the core to remain gravitationally bound. 

Another method to estimate the cloud stability is to calculate the Jeans mass (the 
minimum mass for which a cloud will be bound). The Jeans mass is given by Mj = |7r i?jp 
where Rj is the Jeans length, Rj ~ i'^^) which we take to be the measured radius of the 
cloud, Rdoud- To estimate the mass density, p, we assume a uniform spherical distribution. 
This is clearly not a realistic assumption however the values derived for the number density, 
n, based on this assumption (Table 4) are within a factor of two of the number density 
calculated based on NH3 emission (Table 3). Thus, we assume p probably is good to within 
a factor of 2 or 3. The resulting Jeans mass is given in Table 4. Mdoud/Mj is unity for 
the NH3 gas around G192.16 suggesting that is on the verge of stability while the ratio is 9 
for G192 S3, consistent with the interpretation that the cloud has adequate mass to collapse 
under its own gravitational potential. 
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Table 4: Cloud stability comparison 




G192.16 


G192S3 


Mcloud (Mq) 


0.9 


2.6 


Rdoud (pc) 


~0.07 


0.02 


pt (g cm-3) 


3.7 X 10^20 


3.3 X 10"^^ 


n^ (cni^^) 


8.2 X 10^ 


7.3 X 10^ 


Egrav (ergs) 


8.9 X 10^1 


2.4 X 10^3 


£^t«r6 (ergs) 


1.0 X 10^3 


7.2 X 10^2 


-f^tMrb/^grav 


11.7 


0.3 


Mj (Mo) 


0.9 


0.3 


Mcio«d/Mj 


1 


9 



f Average density assuming a uniform sphere. 



Despite the compact nature of G192S3 and indications that it may be gravitationally 
unstable, no free-free continuum emission from ionized gas has been detected at centimeter 
wavelengths, no H2O maser emission has been detected, and no high- velocity molecular gas 
has been traced to the core which would indicate the presence of a molecular outflow. The 
detection of 850 /im emission toward G192 S3 is consistent with the presence of a warm 
(~ 20 — 30 K), dense molecular core that could be heated by gravitational collapse or an 
embedded protostar. Diffuse 21 /im emission is located up to about 45" south of G192.16 in 
images from the Midcourse Space Experiment (MSX). The emission is at a level of 2 — 3cr 
and there is no strong peak corresponding with G192 S3. The MSX did not detect emission 
near G192S3 at 12 & 14.6 /xm. 

G192S3 has been detected in 1.2 mm continuum emission suggesting that dust is being 
heated by an embedded source (Beuther et al. 2002). Figure 8 presents data from Beuther 
et al. toward G192.16 (05553+1631). Contours have been chosen to illustrate that the 
G192S3 source is detected at a level of Gcr and is clearly separated from G192.16. The 
peak emission at the location of G192S3 is ~ 70 mJy beam^^ in an 11" beam. Following 
the method of Hildebrand (1983), the mass of gas and dust is estimated from the 1.2 mm 



continuum emission using Mgas+dust 



F, D^ 



where D is the distance to the source, F,, 



is the continuum flux density due to thermal dust emission at frequency u, and Bj, is the 
Planck function at temperature T^^. We take k,^ = 0.006( ^^^^^^ )^ cm^ g~^ (Kramer et al. 
1998; Shepherd & Watson 2002). The opacity index P = 1.5 appears to be appropriate 
between wavelengths of 650 microns and 2.7 mm for sub- micron to millimeter-sized grains 
expected in warm molecular clouds and young disks (Pollack et al. 1994). We assume the 
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dust emission is optically thin, the temperature of the dust can be characterized by a single 
value, and the gas-to-dust ratio is 100. Using values of T^ = 30 K and 13 = 1.5, we find 
the mass of gas and dust associated with the 1.2 mm continuum emission peak of G192 S3 
is approximately 4 M©. This value should be considered an upper limit since the G192 S3 
peak fiux density is probably contaminated by the bright, extended emission from the early 
B protostar, G192.16. Never-the-less, this value is within a factor of a few of estimates based 
on NH3 and 850 /xm emission. 

Additional 1.2 mm emission is also detected nearly an arcminute to the SE of the G192.16 
protostar. This is consistent with the presence of 21 /xm and 850 /im emission along a ridge 
to the SE of G192.16, within which G192S3 appears to be embedded. 

We can compare the properties of G192S3 with the well-studied pre-stellar core L1544 
in the Taurus molecular cloud (D = 140 pc) that appears to be on the verge of gravitational 
instability (e.g., Tafalla et al. 1998; Cioleck & Basu 2000; Caselh et al. 2002a, 2002b). 
The L1544 core has a strong 1.3 mm continuum source with a fiux density of 2.3 Jy 
(90 mJy beam~^ peak with a 13" beam or 225 mJy beam~^ with a 22" beam) corresponding 
to a gas+dust mass of 3.2M0 with a density of 1.5 x lO^cm"'^ (Ward-Thompson, Motte,& 
Andre 1999; Caselli et al. 2002a). The N2H"'" line width is ~ 0.3 km s^^. Based on models 
of the chemical evolution in molecular cores, Bergin & Langer (1997) find that NH3 and 
N2H"'" show little depletion in the central, high-density environment of a pre-protostellar 
core because of the low binding energy of nitrogen. Thus, both make excellent tracers of 
very dense regions where depletion is occurring and are expected to have similar line widths. 
Indeed, observations toward several low-mass starless cores in NH3 and N2H"'" show similar 
line widths (~ 0.2 km s~^) for both molecules (Tafalla et al. 2004). Assuming NH3 and 
N2H"'" trace the same regions, one would then expect a similar line width in NH3 as seen in 
N2H+, e.g., ~ 0.3 km s^^ for L1544. 

Although the mass of warm gas and dust in both L1544 and G 192 S3 is similar (a few 
Mq), G192S3 has an NH3 line width that is more than double that expected in L1544 
while the number density is a factor of 10 lower. Massive star formation regions are 
often dynamically active, producing inherently broad line profiles because of large velocity 
dispersions. Despite the increased line width, G192 S3 still appears to be gravitationally 
unstable and may be in the process of collapse. 

Could G192 S3 harbor an L1544-like millimeter continuum peak and not be detected at 
2.7mm? Shepherd et al. (1998) observed the region in which G192 S3 is located at 112 GHz 
with ~ 5" resolution and obtained an RMS of 2.1 mJy beam~^. Scaling the 240 GHz L1544 
emission to 112.3 GHz at a distance of 2kpc {S,y oc u'^/D'^), and assuming that roughly 40% 
of the fiux density would be recovered in the 5" beam, we would expect a peak fiux density 
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of about 1 mJy beam^^ - below the la RMS level in the map. Thus, G192 S3 could contain 
a millimeter continuum peak like that seen in the low- mass core L1544 and not be detected 
at 2.7 mm. On the other hand, the peak flux density in 1.2 mm continuum emission expected 
from L1544 at a distance of 2kpc is 10 mJy beam~^ in an 11" beam. Beuther et al. detect 
a peak flux density of about 70 mJy beam~^ toward G192 S3. Thus, it appears that the 
embedded source in G192S3 is probably somewhat more luminous than L1544. 

This brief comparison illustrates that G192S3 is not likely to harbor a massive 
protostar. The properties are similar to those expected for low-mass pre-stellar or protostellar 
cores. Follow-up, sensitive millimeter continuum observations are necessary to verify this 
interpretation. In addition, N2H+ and HCO"'' observations could provide clues as to the 
evolutionary stage of the cloud collapse (e.g., Lee et al. 2003). 



5. Conclusions 

We have observed the massive star forming region associated with G192. 16-3.84 in 
NH3(1,1), H2O masers, and 1.3 cm continuum emission, and obtained archive data at 850 /xm. 
The dense gas associated with G192.16 is clumpy, optically thin, and has a mass of less than 
IM0. The molecular core appears to be in the process of being dispersed, perhaps by 
the molecular outflow and stellar winds from G192.16. South of G192.16 is a compact, 
optically thick NH3 core (G192S3) with a mass of 2.6 M©. The dense gas is likely to be 
gravitationally bound and may be the next site of star formation in this region. Based on 
the presence of 850 ;um and 1.2 mm continuum emission, G192 S3 may harbor a very young, 
low-mass protostar. 

Water masers trace an ionized jet 0.8" north of G192.16. Masers are also located within 
500 AU of G192.16, their velocity distribution is consistent with the interpretation that 
the maser emission arises in a rotating disk however the current observations do not shed 
signiflcant light on this issue. 
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Figure Captions: 

Figure 1: VLA NH3(1,1) line, 1.3 cm continuum, & H2O maser emission. Top left 
panel shows integrated emission (moment 0) across the main NH3 line from 7.24 to 
4.46 km s~^{vlsr = 5.7 km s^^). NH3 is plotted with thin contours at -2, 2, 3, 4, 5, 
6a (RMS = 2.93 mJy beam~^ km s~^) and in greyscale from la to 17.04 mJy beam~^. 
The 1.3 cm continuum (thick contours) image has an RMS of 0.4 mJy beam~^ and peak 
2.8 mJy beam~^. Contours are plotted at 3, 4, 5a and continue with a spacing of la. The 
synthesized beam of 2.9" x 2.5" at P.A. 44.8° is shown in the bottom right. The top right 
panel zooms in on G192.16; colored symbols represent H2O masers detected by the VLA 
near G192.16 (shown as a filled star symbol). No other H2O masers were detected outside of 
this region. The lower right panel zooms in on a 1.5" region around the massive protostar. 
H2O masers detected by the VLA in this work are shown as small, filled circles with error 
bars. Open circles represent H2O masers detected by Shepherd & Kurtz (1999). Shown 
for reference is a grey torus representing the 1000 AU disk suggested by Shepherd & Kurtz 
based on their H2O maser observations. The small grey ellipse at the center of the torus 
marks the position of the B2 protostar G192.16 and the 130 AU accretion disk proposed by 
SCKOl based on 40 mas resolution observations in 7mm continuum emission. The red and 
blue arrows mark the direction of the large scale molecular outfiow (Shepherd et al. 1998). 
Finally, the black, double arrow north of the torus represents the location of the mixed 
synchrotron and thermal jet detected by Shepherd & Kurtz (1999). Colors of the maser 
symbols represent velocity as shown in the position-velocity plot in the lower left panel. 
The dash-dotted vertical line represents vlsr = 5.7 km s^^. 

Figure 2: The top left panel shows all three epochs of VLB A observations of H2O masers 
(Epoch 1 = triangles. Epoch 2 = squares. Epoch 3 = circles). Top right and bottom left 
panels zoom in on the SW and NE maser clumps. The reference position for which all epochs 
were aligned is shown at postion (0,0). The VLBA masers were aligned to the VLA masers 
as described in the text. Proposed structures in the top left panel are the same as in Fig. 1. 
The bottom right panel presents a position velocity diagram of the masers. 

Figure 3: Left: SCUBA 850/im image of G192.16. The UCHii region associated with the 
B2 protostar is shown as a cross, the southern core, G192 S3, is shown as a plus symbol. 
The center of the image is located at 05^58"'13.87% +16°32'00.1" (J2000). The image RMS 
is ~ 20 mJy beam~^. Contours are plotted as percentages of the peak fiux density of 
2.10 Jy beani"^ (e.g., 10% (lOa) means the contours are 10% the peak fiux density and the 
lowest contour is at lOcr). The beam size of 15" is shown in the lower left of the image. Right: 
An image of G192 S3 after an azimuthally averaged radial intensity profile of G192.16 has 
been subtracted (see text for details). Contours begin at 2a (40 mJy beam^^) and continue 
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at spacings of 20% of the peak flux density (130 mJy beam^^). Compact knots that are 
roughly the size of the beam at the edge of the model-subtracted image are probably due 
to low-level artifacts in the original image. G192 S3 is extended and well within the 120" 
throw, thus, it is likely to be real. Symbols are the same as in the left panel. 

Figure 4: NH3 channel maps. Contours = 2, 3, 4, 5, 6a (RMS = 2.2 mJy beam~^). Velocity 
in km s~^ is indicated in the top left of each panel. The synthesized beam of 2.93" x 2.55" 
at P.A. 44.8° and a scale bar of 0.2 pc are shown in the lower right panel. Top left panel 
shows the zeroth moment map from Fig 1. The filled triangle indicates the position of the 
G192.16 UCH II region. 

Figure 5: NH3(1, 1) spectra toward G192.16 (top) and G192S3(bottom). Main line peak 
= 5.7 km s^^ for G192.16 and 6.7 km s~^ for the southern core. The noise is higher for 
G192S3 because the spectrum was obtained over a smaller area than for G192.16. Thick 
vertical lines below the spectra mark the locations of the main NH3 transition and the first 
two hyperfine components. 

Figure 6: NH3(1, 1) first moment map. Color scale is displayed from dark blue: 4.2 km s~^ 
to dark red: 7.2 km s^^. 

Figure 7: The Spectral Energy Distribution (SED) of G192.16. The filled triangle represents 
the fiux density at 1.3 cm and the filled star is the 850 /im fiux density from this work. Other 
symbols are from: Hughes & MacLeod (1993) - filled square (6 cm); Shepherd & Kurtz (1999) 
- open square (3.6cm); SCKOl - open triangle (7mm); Beuther et al. (2002) - plus symbol; 
Shepherd et al. (1998) - open diamond (3 mm); Williams, Fuller, & Sridharan (2004) - open 
circles (850 & 450 ^um); MSX - filled circles (12.13, 14.65, & 21.3 /im); IRAS - asterisks (100, 
60, 25, & 12yum). IRAS points at 100 and 12 /xm are upper limits. 

Figure 8: 1.2mm continuum emission from Beuther et al. (2002). The image RMS = 
11 mJy beam~^. Contours = 3, 4, 5, 6, 9, 12o" and continue with a spacing of 3a. The 
UCH II region associated with the B2 protostar is shown as a cross, the southern core, 
G192S3, is shown as a plus symbol. 
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